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P_l ■ Abstract 

! We consider self-avoiding walk and percolation in Z^, oriented percolation in Z^xZ. , and 

^ ' the contact process in Z*^, with pD{-) being the coupling function whose range is denoted 

by L < oo. For percolation, for example, each bond {x,y} is occupied with probability 
pD{y — x). The above models are known to exhibit a phase transition when the parameter p 
CN \ varies around a model-dependent critical point pc- We investigate the value of when d > 6 

for percolation and d > A for the other models, and L ^ 1. We prove in a unified way that 
1/-^ , Pc = 1 + C{D) + 0{L^'^'^), where the universal term 1 is the mean- field critical value, and 

I the model-dependent term C{D) = 0{L~'^) is written explicitly in terms of the random walk 

■ transition probability D. We also use this result to prove that Pc = 1 + cL~'^ + 0{L~'^~^), 

Tj- . where c is a model-dependent constant. Our proof is based on the lace expansion for each of 

' these models. 



1 Introduction and main results 

Self-avoiding walk, percolation, and the contact process are well-known models that exhibit critical 
phenomena. For percolation in two or higher dimensions, for example, there exists a percolation 
threshold pj!'' such that there is almost surely no infinite cluster for p < p^", while for p > pl" there 
is almost surely a unique infinite cluster. As p j p^", the average cluster size and the correlation 
length diverge. The precise value of pl" depends on the details of the model, and is only explicitly 
known in a few cases, such as for two-dimensional nearest-neighbor bond percolation . 

In this paper, we will consider self-avoiding walk, percolation, oriented percolation and the 
contact process, where the interaction range L is taken to be large. When L ^ 1, the interaction 
in the considered models is relatively weak, and therefore the critical values can be expected to be 
close to the critical value 1 of the respective mean- field models, i.e., random walk and branching 
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random walk. We study the difference of the critical values and 1 for the above four models as 
L — >• oo. It turns out that, above the respective upper critical dimensions, we can write this 
difference to leading order as simple functions of the underlying random walk. 

1.1 Models 

First, we define the models. A self-avoiding walk is a path uj in the dimensional integer lattice 
1/ with uj{i) 7^ uj{i) for every distinct j G {0, 1, ... , \uj\}. We also take the zero-step walk into 
account. We define the weight of a non-zero path ui by 

D{uo{t)-uo{t-l)), (1.1) 

i=l 

where D is a probability distribution on I/, and let Wp{uj) = 1 if |ci;| = 0. We suppose that D 
is symmetric with respect to the lattice symmetries and that D{o) = 0, where o is the origin in 
Z"*. A more detailed definition will be given below. The self- avoiding walk two-point function is 
defined by 

r-(x)= ^pH' (1-2) 

W.O rX 

saw 

where the sum is over all self-avoiding paths from o to x. It is known (see, e.g., [22]) that there is 
a critical value pl"" such that 

x7 = Yl ^r(^) (1-3) 

is finite if and only if p < pl"" and diverges as p t pl^. 

For percolation, each bond {x,y} is occupied with probability pZ}(y — x) and vacant with prob- 
ability l—pD{y — x), independently of the other bonds, where p G [0, ||-D||oo^]. Since D{x) = 1, 
the percolation parameter p is the expected number of occupied bonds per site. We denote by 
Pp the probability distribution of the bond variables. We say that x is connected to y, and write 
X < — > y, if either x = y or there is a path of occupied bonds between x and y. The percolation 
two-point function and its sum over Z'^ are denoted by 

r-(x) = Pp(o ^ x), xr = E ^r(^)- (1-4) 

Similarly to self-avoiding walk, there is a critical value p^" such that Xp" is finite if and only if 
p < pl" and diverges as p t pl° (see, e.g., |5|). 

Oriented percolation is a time-directed version of percolation. Each bond ((x,t), {y,t + 1)) is 
an ordered pair of sites in Z'^ x Z+, and is occupied with probability p D{y — x) and vacant with 
probability 1 —pD{y — x), independently of the other bonds, where p G [0, We say that 

(x, s) is connected to (y, t), and write (x, s) — > {y, t), if either (x, s) = {y, t) or there is an oriented 



WpiLu)=p^' 
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path of occupied bonds from (x, s) to {y,t). Let Pp be the probabihty distribution of the bond 
variables. The oriented percolation two-point function and its sum over Z'' x Z+ are denoted by 

T;-{x,t)=F,{{o,o)-^{x,t)), x°;=Yl (1-5) 

Also oriented percolation exhibits a phase transition such that < oo if and only if p is less than 
the critical value and that x°p | cxd as p t (see, e.g., PJ. 

The contact process is a model of the spread of an infection in Z*^, and is a continuous-time 
version of oriented percolation in Z'' x ]R_,_. We now describe a graphical representation for the 
contact process. Along each time line {x} x ]R_|_, where x G Z^, we place points according to a 
Poisson process with intensity 1, independently of the other time lines. For each ordered pair 
of distinct time lines from {x} x ]R_|_ to {y} x ]R_|_, we place oriented bonds {{x,t), {y,t)), t > 0, 
according to a Poisson process with intensity pD{y — x), independently of the other Poisson 
processes, where the parameter p > is the infection rate. We say that (x, s) is connected to 
(?/,t), and write (x, s) — > (l/, if either (x, s) = (|/,t) or there is an oriented path in 1/ x ]R_,_ 
from (x, s) to (y, t) using the Poisson bonds and time-line segments traversed in the increasing- 
time direction without traversing the Poisson points. Let Pp be the corresponding probability 
distribution. We denote the contact process two-point function and its integro-sum over 1/ x ]R_|_ 
by 

/■oo 

r7(x,t)=Pp((o,0)^(x,t)), x7= dtJ2Vi^^t)- (1-6) 

Jo 



Again there is a critical value p"^ such that x^p is finite if and only ii p < p"^ and diverges as p t P^c 
(see, e.g., EH)- 

We will omit the superscript referring to the precise model, and write pc when referring to the 
critical values in all models simultaneously. The goal in this paper is to study Pc when the range 
L of D is sufficiently large. In the proofs, we will have versions of D in mind which are such that 
L^D{Lx) is a discrete approximation of a function on M'^. We will formalize this assumption on 
D in the following definition: 

Definition 1.1. Let hhe a. probability distribution over ]R^\{o}, which is invariant under rotations 
by 7r/2 and refiections in the coordination hyperplanes. We suppose that h is piecewise continuous, 
so that Jjgd h{x) dfix = 1 can be approximated by the Riemann sum L^'^ J2xeZ'' h{x/L). Then, we 
define 

hjx/L) 

We will make heavy use of results proved elsewhere for the models under consideration. For 
these results, some further assumptions are made on of which we now list the most important 
ones. We require that there exist c > 0, C < oo, 77 G (0, 1) such that 

sup D(x) < CL^^ r] A{cL'^\k\'^) <1- D{k) <2-r], (L8) 
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where D{k) = ^^^^d D{x) e^^'^ and \k\'^ = Yl'j=i^'j- There are a few more minor requirements 
that depend on the precise model under investigation. For details, see [TT] for percolation and 
ITH| ITT)| IT7j for the other three models, for which the requirements are virtually identical. A 
simple example of D, where all the above assumptions are satisfied, is 

M0<M^<L} 

for which h{x) = 2~'^ if < ||a;||oo < 1 and h{x) = otherwise. 

We denote hj D * G the convolution of D and a function G in Z'^, and by D*" the n-fold 
convolution of D in Z"*, where we define D*^{x) = 6o^x- We will frequently use 

D*-ix)<6o,Jo,x + ^^^, (1.10) 

where 

13 = L-'^. (1.11) 
The inequality p.l()|l is a consequence of ()1.8j) . as we will show in Appendix 1X1 



1.2 Main results 

Let dc denote the respective upper critical dimensions, i.e., dc = Q for percolation and dc = ^ for 
the other three models. In this paper, we investigate the respective critical values when d > dc 
and L ^ 1, in a unified fashion. 

Theorem 1.1. For each model with d > dc, as L oo, 

oo 

p-, = 1 + 5^ D™(o) + 0(/?2), (1.12) 

n=2 
^ oo 

pI^ = 1 + -J2d*'^o) + 0{(3'), (1.13) 

n=2 

oo 

= 1 + D*\o) + - Y^{n + 1) D™(o) + 0{f3^). (1.14) 

n=3 

The universal term 1 is the critical value for the mean-field models (random walk and branching 
random walk). Note that, by p.lOj) . the model-dependent terms in fll.l2|) - ()1.14j) are 0(/3). In 
Section 11.31 we will intuitively explain why the model-dependent terms have the above respective 
forms. 

We next compute the dependence on f3 more explicitly, and compute the coefficients of /3 in 
Pc — 1 explicitly. For this, we let U be the uniform probability distribution over [—1, 1]*^ C M"^, i.e., 
for X e M'^, 

f/(x) = 2~"'l{||x||^<i}, (1.15) 

and denote by U*"' the n-fold convolution of U in M'^. Then, the leading order coefficient in /3 for 
Pc is given in the following theorem: 
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Theorem 1.2. Fix D as in ()1.9p . and let d> d^. As L oo, 



oo 



Pc, Pc 



l + ;5^f/-(o)+0(/3L-i), 



(1.16) 



n=2 



l + |^f/-2n^0) + 0(/?L-l) 



(1.17) 



n=2 



1 + /5 f/*2(^)_^_^(^_^^)^.n(^) +0(/3L-l). 



(1.18) 



n=3 



We now comment on the relation between the asymptotics in Theorems II . iWT!^ The advantage 
of Theorem 11.21 is that it is more concrete, and the continuum hmit of the critical points appears 
explicitly. However, the error term in Theorem 11.11 is 0(/3^), while in Theorem 11.21 it is equal to 
0{(3L~^), which is much larger. In order to compute the critical value more precisely, Theorem ll.il 
gives a much more powerful result, at the expense of having to compute the random walk terms 
appearing in its statement. In principle, it should be possible to compute the coefficients of 
PL~^, PL'"^, . . . , PL~'^~^^, but this requires a substantial amount of work. Finally, it should be 
possible to compute the random walk sums in Theorem 11.11 for other examples than the one in 
()1.9|) . but we refrain from doing so. 

We now summarize previous results on the critical values. We start with self-avoiding walk. 
Penrose's result in [25j implies that the critical value for self-avoiding walk defined by ()1.9j) with 
L ^ 1 satisfies 



for some /3-independent constants c, c', c". For spread-out lattice trees, a related result with a 
different leading term, namely = e~^, was also obtained in [22] • For d > 4, Madras and Slade [23 
Corollary 6.2.7] improved ()1.19|) to pl"" = 1 + 0{P). In [I3[T7j, this result was extended to more 
general D as defined in Definition ll.il We will rely on the results in jTHlEI, whose proof is based 
on the lace expansion and a generalized inductive approach. We will also use the lace expansion 
to derive the expression of the 0(/3) term in p.l2|) . 

For percolation, the best previous result is pl" = 1 + 0(/5^/'^~'^) for d > 6 and L ^ 1, where 
e > is an arbitrarily small number |10j. However, if we combine Lemma l3. II proved below and 
the estimates for the lace expansion in fl], then we obtain the better estimate p^" = 1 + 0(/3). 
The result in p.l4p . which is also obtained by an application of the lace expansion, identifies the 
expression of this 0(/3) term. 

When d > 4 and L ^ 1, both and pl^ were proved to be 1 -|- 0{f3) ^2 ^1 CSl- Similarly 
to self-avoiding walk, the proofs of these results rely on the lace expansion and an adaptation of 
the inductive approach. The contact process defined in terms of D of p.9j) was first considered by 



{1, a d> 3, 

l + c'(3\og(3-\ if d = 2, 
l + c"/3^/^ if rf= 1, 



(1.19) 
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Bramson, Durrett and Swindle and they proved that, as L — > oo, 

if c/ > 3, 

p--lx/(/3)^<(/31ogr\ ^id = 2, (1.20) 

if = 1, 

where p^"" — 1 x /(/3) means that the ratio {pl^ — 1)/ f{(3) is bounded away from zero and infinity. 
Later, Durrett and Perkins [7j proved that 

Our result p.lfij) in Theorem 11.21 is stronger when ci > 4 in the sense that not only the coefficient 
of /3, but also the speed of convergence in p.2ip is identified. In [T3j, we also obtained certain lace 
expansion results for a local mean-field limit, where the range and time grow large simultaneously, 
for the contact process in < 4, and we expect that these results could be used to prove a stronger 
version of p. 211) for = 3, 4, as well as for oriented percolation. However, this will need serious 
work using block constructions as used in [Tj. 

We expect that ()1.12|) - ()1.14p remain valid ioi d = — 1 and dc when we change 0(/3^) to 
o(/3). As mentioned above, this is the case for the contact process jZ]. When d < dc — 2, the 
second terms in p.l2j) - ()1.14j) diverge, so that Theorem 11.11 cannot hold. However, we expect that 
the asymptotics of the critical point will, as for the contact process, again be described by the 
divergence of the sums in p.l2|) - ()1.14|) . 

When d > dc, we expect that the 0(/5^) terms could be identified in terms of D as well, using 
a similar method as in this paper, but to do so will require a serious amount of work. 

A related problem is to obtain the asymptotics of the critical points for the nearest-neighbor 
models, when D{x) = {2d)~'^l{\x\=i} and d —>■ oo. In [12j, pl"" was proved to have an asymptotic 
expansion into powers of {2d)~^, and the first six coefficients were obtained. For unoriented perco- 
lation, the first three coefficients were computed in [T2j and [IH], but the proof of the asymptotic 
expansion only appeared in JHl- The proofs of these results are again based on the lace expansion. 
For nearest-neighbour oriented percolation and the nearest-neighbour contact process, it is proved 
that = 1 + 0{d~^) (see 0) and pl^ = 1 + 0{d~^) (see, e.g., |22]), using different methods. 



1.3 Overview of the proof 

To prove Theorem I l.H we will apply the lace expansion (see, e.g., [TTl[l3[[T21|2ni|2lj). For example, 
the lace expansion for self-avoiding walk gives the recurrence relation 

Tpix) = 6o,. + 5^[p/^(t;) + H-(^;)] r-(a; - v), (1.22) 

V 

where Il''p{x) is a certain expansion coefficient. It was proved in [13 CZj that 11°'' = J2x^7^^^ ~ 
0{j3) for p < pl"", if c? > 4 and L ^ 1 (see Section 12)). Summing both sides of p.22j) over x G Z'' 
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and solving the resulting equation in terms of x'pi we obtain 

X- = (l-p-fl-)^\ (1.23) 

and thus 

Pr = l-n;|a. (1.24) 

To estimate p'^, we thus need to investigate Ilpsa. We will prove that, since p'^ = 1 + 0(/3), we can 
replace flpsa by fl^^^ up to an error of order 0(/3^). When p = 1, the exponentially growing factor 
pl*^! in (ll.lj) does not play any role, and Ilf can be investigated in terms of random walks. This is 
the key ingredient for the proof of Theorem 11.11 

The strategy for percolation models is the same as above. There is a similar recursion relation to 
()1.22j) . with some mo del- dependent expansion coefficient 11^(0;). Therefore, to obtain the formulae 
in Theorem I l.H we will have to investigate Hi = Ili{x), again in terms of random walks. 

As we will explain in Sections |2H3 ftf and Hi can be described by an alternating sum of a 
model-dependent sequence 7i['^^ for > 0, where fr['^^ for > 1 decays as for all models. 
For self-avoiding walk, frf^ equals zero, while vrj"' for percolation models is nearly a half of Trj^''. 
(This is why we have the factor ^ in (jl.l3j) - (jl.l4j) .) Therefore, roughly speaking, we only need to 
investigate ttJ^^ to obtain p.l2|) - ()1.14|) . We will show later that the diagrammatic interpretation 
of Tii^ for self-avoiding walk is a single random walk taking more than one step and going back to 
the starting point (cf., p.l2|l ). while the diagrammatic interpretation of tt^i^ for percolation models 
is that two random walks, at least one of which is non-vanishing, meet at some point. Therefore, 
the correction to the mean-field value 1 are related to random walk loops. 

For loops in the time-oriented percolation models, the lengths in the time-increasing direction of 
these two walks have to be equal (which explains the sum over even convolution powers in ()1.13|l ). 
while for unoriented percolation this is not the case (which explain the sum over all powers and 
the factor n + 1 in ()1.14|l ). 

For the contact process, the two paths are continuous time random walk paths, for which the 
number of convolution powers of D is equal to the number of spatial steps made by the random 
walk up to a given time, which has a Poisson distribution. Therefore, the sum over the convolution 
powers of D is not restricted to even powers, and we see that the correction to the mean-field value 
for the contact process and oriented percolation are different. For the contact process, it will turn 
out that also the factor 1 in ()1.13|) disappears, which is due to the fact that the two walks are 
in fact avoiding each other, and which will be explained in more detail in Section 13.11 This is an 
intuitive explanation of the mo del- dependent terms in (|1.12j) - (jl.l4j) . 

We organize the rest of this paper as follows. We begin with self-avoiding walk in Section |21 
and explain the key steps to estimate pl^. Following the same steps, we estimate p"^ and pl^ in 
Section 13.11 and p^" in Section 13.21 Finally, we prove an extension of ()1.10|) in Appendix |X1 and 
Theorem 11.21 in Appendix iBl 
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2 Critical point for self-avoiding walk 

In this section, we prove (jl.l2j) . using (jl.24|) . Throughout this section, we will omit the superscript 
"sa" and write, e.g., Pc = pl'' and Tip = 11^. 

Before computing the asymptotics of IIp^ in ()1.24|1 . we first note that Pc > 1- This is because 

the removal of the self-avoidance constraint in p.2|) results in J2u)-o »x^p('^)' whose sum over 

X ^li^ equals (1 — for any p <\. For self-avoiding walk, 

oo 

np(x) = (2.1) 

where, e.g., vr^^'(x) is a "1-loop diagram" at the origin [23]: 

T^'p\x) = 5o,. {pD * Tp){o) = 6o,. Yl ^p(^) ^(^)' (2.2) 

LU'.O >0 

\lo\>1 

where I{uj) = 1 if there are no self-intersection points except for u;(0) = cjdu;!), otherwise I{uj) = 0. 
For d > 4 and L ^ 1, it was proved in p^J that, for tt^'^' = J2x ^p'^'(^)) have 

frir'<0{(3f, pdpllp<0{(3), (2.3) 

for allp < Pc and > 1. Together with ()1.24|) and ()2.H) . we immediately obtain that Pc = l+0(/5). 
Moreover, by the mean- value theorem, there is a p G (l,Pc) such that 

= 1 - fli - (flp, -U,) = l-U,- (p, - 1) dpflp = 1 + 7r« + 0{(3'), (2.4) 

where 

oo 

^':'= E W,{uj)I{uj) = J2d*''{o)- E W,{u;)[1-I{u)]. (2.5) 

uj:o >o 71=9 ^'-O >o 

\uj\>l \uj\>l 

To complete the proof of ()1.12j) . it thus suffices to prove that the second term in the right-hand 
side of (123) is 0{p'^) if > 4. We first note that /( uj) is an indicator function. If I{uj) = 0, so 
that 1 — I{uj) = 1, then there must be a pair {s,t} ^ {0, \ijjW with < s < t < |c<j| such that 
ioi^s) = uj{t). Denoting the parts of uj corresponding to these three time intervals by tUj, z = 1, 2, 3, 
respectively, we obtain 

E W^iM [1-/M]< E E E flWM) = iD*G*')io){D*'*G){o), (2.6) 



uj:o >o 



g^d uJi,LJ[i:o tx UI2-X >x i=l 



I^I>1 |<^i| + |lJ3|>1 |w2|>1 

where G{x) = T.7=o and {D*^ * G){o) = Er=2^*"(o) = 0{(3) ii d > 2. Moreover, by 

(mni), 

oo 

{D * G*2)(o) = E ^^™(o) = 0{(3) (2.7) 

n=l 

if d > 4. This completes the proof of p.l2|) for self-avoiding walk. □ 



8 



3 Critical points for percolation models 



In this section, we compute the asymptotics of the critical values for the other three models, and 
thus complete the proof of Theorem ll.ll 

To discuss oriented percolation and the contact process simultaneously, it is convenient to 
introduce the following oriented percolation on Z*^ x which is the time-discretized contact 
process with a discretization parameter e e (0, 1]. A bond is a directed pair {y.t + e)) of 

sites in 1/ x £:Z+. Each bond is either occupied or vacant, independently of the other bonds, and 
a bond ((x, t), (?/, t + e)) is occupied with probability 

Qpiy-x) = \ jf/ . (3.1) 

provided that sup^gp(a;) < 1. In this notation, the model with £ = 1 is the usual oriented 
percolation model as defined in Section [HTJ and the weak limit as e | is the contact process 
Similarly to oriented percolation with e = 1, for each e G (0, 1], there is a critical value p^"^ for 
every e G (0, 1], such that p^}^ = and hmej^oPc^'' = P^c EEj- We will call the model with e G (0, 1] 
the time-discretized contact process. 

To summarise notation for percolation and the time-discretized contact process, we will write 
K = !/• for percolation and A = Z*^ x eZ+ for oriented percolation. For notational convenience, 
we will take e = 1 for percolation. We will also use bold letters to represent elements of A. For 
example, o = o, a; = a; for percolation, and o = (o, 0), a; = {x,t) for the time-discretized contact 
process. For a bond h = {u, v), we write h = u and b = v. We also omit the superscripts e, pe, op 
and cp, if no confusion can arise. 

As mentioned in Section ^ the lace expansion for percolation models takes a similar form as 
in (fT^ . and reads (see, e.g., [TTllTlj l 

Tp{x) = [5a,x + np(x)] + ^ [6o,u + ^p{u)] qp{v - u) Tp{x - v). (3.2) 

In particular, qp{v — u)= pD{v — u) for percolation and oriented percolation for which e = 1. (To 
unify notation, we recall that we regard unoriented percolation as a model with e = 1.) The lace 
expansion coefficient Hp (a?) equals 

oo 

Ilp{x) = 5^(-l)^<'(^), (3.3) 

Af=0 

where 7r^'^^{x), N > 0, are model-dependent diagram functions. The result of the lace expansion 
will be explained in Sections I3.1H3. 21 For the time-discretized contact process with e E {0,1], d > dc 
and L 3> 1, it has been proved [El CHj that lip = eJ^xeA^pi^) 0{/3)e'^ for all p < pc- The 
same estimate is proved to hold for unoriented percolation (with e = 1), using the lace expansion 
in jll] and Lemma f3.1l proved below in Section l3?2l 
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As in the derivation of ()1.23|) . solving ()3.2|) in terms of Xp = £ J2xeA '^p(^) gives 

V = L±i55 (3 4) 

and thus, equating the denominator to zero, 

Pc = l-^K-{Pc-l)l%.. (3.5) 

This expression holds uniformly in e. We will use it to compute pl'^ and pl" by taking e = 1 and p^F 
by taking the limit when e I [221; respectively. In particular, the third term is 0(/3^) when e = 1, 
and it has no contribution in the limit e [ 0. Therefore, we are left to prove that, apart from an 
error term of order 0(/5^), the second term in (j3.5|) equals the second term in (jl.l2p when e i 0, 
and equals the second term in ()1.13|) for oriented percolation and that in ()1.14p for (unoriented) 
percolation when e = 1. We again note that p^"' > 1, since Xp < e X]^o Ip^i^) = (1 ^ P)~^ 
p < 1. In addition, similarly to (|l.iup . when p = 1 and e < 1, we have 

g-(x) ^ gr(x) < (1 - sT + J^^^r (3-6) 

We will prove ()3.6|1 in Appendix El Note that when e = 1, ()3.6|1 reduces to ()1.10|) . 

To complete the proof of Theorem II. H we investigate IIp^ for oriented percolation and the 
contact process in Section 13.11 and for unoriented percolation in Section | 



3.1 Asymptotics of pl^ and p"^ 

In this section, we investigate IIp^ for the discretized contact process, and derive ()1.13|) for oriented 
percolation (i.e., s = 1) and p.l2|) for the contact process (i.e., e I 0). 

To describe the diagram functions '7r^^'(a;), > 0, we need some definitions. We say that x 
is doubly connected to y, if either x = y or there are at least two nonzero bond-disjoint occupied 
paths from x to y. Following the notation in as closely as possible, we denote this event by 
X =^ y, and define 

^p ' = ^ E = ^x)- 5^,^. (3.7) 

xeA 

If o is connected but not doubly connected to x, there is a pivotal bond b = [b, b) for o — > x such 
that both o — > b and b — > x occur, and that o — > x occurs if and only if b is set occupied. For 
A C A, we say that y is connected to x through A when every occupied path from y — > x has 
at least one bond with an endpoint in A. We define E{b, x; A) to be the event that b is occupied, 
that b — > X through A, and that there are no pivotal bonds b' for b — > x such that b — > 6' 
through A. Let C^{o) be the set of vertices in A connected from o without using b. Then, 

7r« =£^<H^), where 7r^;\x) = ^Pp(o^ 6; Eib,x;C\o))). (3.8) 

xeA b 
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The higher order diagram functions TT^p^^{x), N >2, are defined in a similar way, but are irrelevant 
in this paper (see ^1 Section 3] for a complete definition, with slightly different notation). 
For d > 4 and L 1, it was proved in [14^ that, for tt^^' = e J2xeA '^p^K^)^ have 

TT^"' < OiPf^'e', pd,fl, < 0{P)e', (3.9) 

for all p < Pc and > 0. Together with ()3.3p and ()3.5|) . we obtain pc = 1 + 0(/5). Moreover, by 
the mean- value theorem, there is a p G (l,Pc) such that 



= 1 - - (Pc - 1)-% = 1 - ^Hi - (pe - l)^dpU, + 0{[3') e 

1 1 



+ ^7rf^+0(/52). (3.10) 

To prove (jl.l2|) - (|1.13p . it thus suffices to investigate irf and tx'^K 
Analysis of 7if\ We prove 

1^,0) f=^Er=2^*'"(o) + 0(/32), for 5 = 1, 

' l^E"=2^*"(o) + C'(/32), when ejO. ^' ' 

Recall ()3.7|) . To describe a double connection by a pair of two random walk paths, we order 
the support of D in an arbitrary but fixed manner. For x, y in the support of D, we write x -< y if 
X is lower than y in that order. For a pair of paths consisting of bonds in A, a; = (6i, . . . , b^f) and 
u' = {b[, . . . , 6^) with = and b^, = b'^, we say that uj is lower than tu', denoted by a; -< u', if 
at the first time n G {1, . . . , A^} when uj is incompatible with a;' (therefore 6j = 6^ for all i < n) we 
have bn -< We also say that uj2 is higher than cui. 

A path LJ = (6i, . . . , 6|^|) is said to be occupied if all bonds along u are occupied. We define 
E^{uj) to be the event that u is the lowest occupied path among all occupied paths from b^ to b\;^\, 
and that there is another occupied path u' from b^ to b\i^\ which is bond-disjoint from u (denoted 
hj UJ n uj' = 0) . Given a path uj, we also define E^{uj'] uj) to be the event that u' is the highest 
occupied path among all occupied paths from b^ to that are bond-disjoint from uj. Such an 
occupied path uj' exists on {b^ =^ ^ -^-<(^) by definition. 

Using the above notation, we have, for x ^ o, 

{o x} = \^ [u!i,u!2 occupied; E^(u!i) n Ey(uj2;uJi)} ■ (3-12) 

UJl,UJ2'.0 >x 

We define the right-hand side to be empty if x = o. Then, 

Pi (cui, co'2 occupied; -E'^(to'i) n -E^(ciJ2; cc'i)) . (3.13) 

X£A ^1,^2-0 >X 

LUinu}2 = l2> 
iJ\^LU2 
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Since Pi is a product measure, if we ignore E^{ui) fl Ey{u!2; uji), then we obtain 
Pi(c<Ji, ct;2 occupied) 



UJl,UJ2-0 >x 



= q{u) q{v) q{x — y) q{x — z) Pi (cji occupied) Pi (0^2 occupied), (3-14) 

u,v:u-iv u)\:u >y 

where u = {u,e), v = {v,e), y = {y,t — e), z = {z,t — e), and q{x) = qi{x) (cf., ()3.6|) ). By an 
inclusion-exclusion relation, the correction is bounded by 

[Pi (cui , ci;2 occupied; E^^uiY) + Pi (cui , ^2 occupied; £'^(072; (x^i)'^)] . 

uJi,iU2:o 

a;inw2=0 

We will prove below that, for E equal to E^{uJi) or Ey{uj2; uji), 

eJ2 PiK, ^2 occupied; E=) = 0(/52) (3.I5) 

Lu-ir\ijj2=0 

We investigate ()3.14|) to obtain the expression of 0{(3) from ()3.13p . If we ignore the restriction 
tui n c<j2 = 0, then we obtain 

q{u) q{v) q{x - y) q{x - z) q*^'^'-^\y - u) q<'''-^\z - v), (3.16) 

where tje G [2, 00) fl Z_|_. We will prove below that the correction satisfies 
e Y^ li'^) ~ y) ~ ^) Pi(t^i occupied) Pi (0^2 occupied) = O{0^) e^. (3.17) 

£cgA u,v:u^v ujiiu >y 

1^10(^27^0 

Therefore, we only need to consider the contribution to ()3.13|) from ()3.16p . By changing variables 
as y' = X — y and z' = x — z and using the symmetry between u ~< v and u >~ v, the sum of ()3.16p 
over X G Z*^ equals 



J2 Qiu) q{v) q{y') q{z') J] g*(*/--2) {x - y' - u) g*(*/--2) {x - z' - v) 

X 

\ Y q{u)q{v)q{y)q{z)q<'"^~'\v + z-y-u). (3.18) 



y',z':y'^z' 



y,z-y¥'z 



Recall 1)3.11) . Since there is at most one temporal (or vertical) bond growing out of every site in 
A, we must have q{u) = eD{u) or q{v) = 6D{v), so that we obtain at least one factor of e. By the 
same reason, we should have q{y) = eD{y) or q{z) = eD{z), so that we obtain a second factor of 
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e. Therefore, the number of combinations for the product of four factors of q in ()3.18p is nine: one 
combination is proportional to e^, four others are proportional to (1 — e)e^, and the remaining 
four are proportional to (1 — e^. Only the first case arises for oriented percolation for which 
£ = 1 , while only the third case arises for the contact process for which £ i 0, respectively. 

We first complete the proof of (j3.1ip for oriented percolation. When e = 1, and using inclusion- 
exclusion on the restrictions u ^ v and y ^ the sum of ()3.18|) over t>2 equals 

^ oo _ oo 

- Yl Diu)Div)Diy)Diz)J2D*^''''\v + z-y-u) + 0{f3') = -J2D*'\o) + 0{(3'), 

u,v,y,z t=2 t=2 

(3.19) 

where we use fjl.8|) to obtain an error of order 0(/3^) that comes from contributions where u = v 
OT y = z. 

For the contact process, for which e I 0, the leading contribution is due to the four combinations 
of order (1 — e)"^ e"^ mentioned above, where either m or f is o, and either y or z is o. Therefore, 
the coefficient of (1 - e)^^^ in (nTTHll is 
1 



J2d{u) D{y) q*^'"'-'\-y - u) + Y.D{u) D{z) q<^"^-'\z - u) 



v,y 



2(D'''tq'^''l'-''>)(o). 



Summing this expression (multiplied by e) over t/e E [2, oo) fi Z+ gives 

d'^k 2b{kf 

[1 — fc -t- biy\Kj]^ = j 

^ f A,^^ = yD-{o), (3.20) 

where 0,^ = 7r,7r]'^. This completes the proof of ()3.1H) . □ 

Analysis of tx^^\ We prove that \tIi^ is asymptotically twice as large as the right-hand side of 
dSHD: 

l^a, f=Er=2^*'"(«)+0(/?^), for 6 = 1, 

' 1^2Er=2^™(«)+0(/5'), when ejO. ^ ' ' 

For a bond 6, let {h => x} be the event that b is occupied and b x. We define {u — > b} 
and a joint event {u — > b =^ x} similarly. For events Ei and E2, we denote by Ei o E2 the event 
that El and E2 occur disjointly, i.e., using disjoint bond sets of bonds (see e.g., [HI Section 2.3]). 
Recalling ()3.8|) and distinguishing between b = o and 6 7^ o, we can rewrite tt^^^ as 

7rf^=£ Yl Pi({(o,^x)^x}o{o^x})+£5^ ^ Pi (0^6; E{b,x;C\o))) 

-6 J2 ^i{{{io,u) ^ x} o {o ^ x}} \ Eiio,u),x;C^''^-\o))y (3.22) 
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We will extract the leading contribution from the first term. Note that {(o, u) — > x} o {o — > x} 
is almost identical to {o =^ x} = {o — > x} o {o — > x}. However, the symmetry between 
the two connections from o to a; is lost in the former event, due to the bond (o, u). We will use 
this symmetry breaking in a convenient manner. Recall that below ()3.1H) . the support of D was 
ordered in an arbitrary way. Now, instead, we choose the ordering such that, for u = {u,6), the 
element u in the support of D is minimal. This will ensure that the lowest occupied path ui from 
o to a; will use the bond {o,u). We also write and E'!^{lj2]^^i) for E^{uJi) and Ey{LJ2]uJi) 

in this n-dependent ordering. Therefore, (cf., ()H.12j) ). 

{{o,u) — > x} o {o — > x} = [J {uJi,iU2 occupied; E'!^{uJi) n E'^{iU2',^i)}, (3.23) 

uji:(o,u) >£C 

0^2:0 — >cc 

ujinuj2=0 

and its contribution to ()3.22|) is 

^ Pi (cui, W2 occupied; ^"(t^i) n E^{uj2] uji)) , (3.24) 

tcgA (x'i;(o,it) >£C 

LU2'.0 >X 

where uoi : (o, u) — > a; is a path from o to a; starting by the bond (o, u). Ignoring the condition 
E'^{uji) n E^{uj2;uji) as in ()3.14j) and following the same strategy as in estimating 7tf\ we obtain 
the main contribution to 1)3.211) . The leading term of ^vrj^' is twice as large as that of 
because the symmetry is broken and we do not obtain the factor | as in ()3.18j) (cf., ()3.13j) and 

(EH). 

To complete the proof of ()3.21)1 . it suffices to show that the second and third terms in ()3.22j) 
are both 0{(3'^) e^. The event in the second term of ()3.22|) implies the existence of y & A such that 
{o — > y — > h} o {o — > 6} and {y — > x} o {h — > x} occur disjointly. Let uji denote a path 
from o to a; through y, UJ2 denote another path from o to a; via the the bond h with h = z, and 
uj^ denote another path from y io z. Then, the second term in ()3.22p is bounded by 

3 

x,v,zeA ^i-o — *y — 1=1 

z^ox '^2:0 >z *x 

since Pi is a product measure. The third term in ()3.22|) is also bounded by the above expression. 
This is because the event in the third term in ()3.22)1 implies existence of ?/ G A and a pivotal bond 
h = (z, ■ ) for u — > X such that {o — > y — > as}, {(o, u) — > h — > x} and {y — > z} occur 
disjointly. We thus obtain ()3.25j) by the same random walk representation. 

Therefore, it is sufficient to prove that ()3.25|) is bounded by 0(/3^) e^. When e = 1, we simply 
ignore the restriction uJi H uj = 0, i j , and apply the Gaussian bound ()1.10|1 to the part of ui 
from y to X and to the part of UJ2 from o to z. Since y ^ x and z ^ o, the term So,x in (jl.lOj) 
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does not contribute, so that ()3.25|) is bounded by 



0<s<s'<t 

where s,s' are the time variables of y and z, respectively. When e < 1, we use the restriction 
uJinujj = 0, i ^ j, to extract factors of q with pairwise different arguments, as in ()H.14j) . out of the 
four intersection points o, y, z and x. As explained above ()3.19|) . each pair gives rise to a factor 
e, and we obtain a total factor e"^. With the help of ()3.6|) . ()3.25|1 with e < 1 is bounded by 
times the left-hand side of fl3.2fi|) with the region of summation being replaced by eZ+. This is 
further bounded by e^, since the sum over t, s, s' G eZ^ eats up a factor for the Riemann 

sum approximation. This completes the proof of ()3.21|1 . □ 

Proof of (I3.15j) . We only consider the case = E'^{uJiY, which is the event that there is an 
7] -< Ui from o to x, which must share at least one step with uJi, such that E'^ij]) occurs; the other 
case E = Ey{uj2] <^i) can be estimated in a similar way. Let be the part of r] from the point, say 
y, where r] starts disagreeing from uji until it hits ui or UJ2 at z. Since Pi is a product measure, 
()3.15p is bounded by 

3 

e ^ ^ (llzet.'Aiy}} + l{2Gcj2}) JJ^ Pi (coij occupied). 
x,v,zeA ^i-o — *y — *^ 1=1 

z^ox '^2:0 >x 

U>inLJj=0, iy^j 

Since the contribution from l{zeuj2} is equal to ()3.25|) . we only need to investigate the contribu- 
tion due to the other indicator l{ze<^i\{y}}. We again discuss the case e = 1 first, and then adapt 
the argument to the case e < 1, as done below ()3.26|) . When e = 1, we ignore the restriction 
Ui n Uj = 0, i j , and apply (jl.lOj) to the probability of UJ2 and being occupied. By denoting 
the time variables of y and z by s and s' respectively, the contribution from l{zeuji\{y}} is bounded 
by 

0<s<s'<t 

When e < 1, we use the restriction tUj ft cj^ = 0, i 7^ j, along each of the four intersection points 
and obtain the eight factors of q with pairwise different arguments. Following the argument below 
()3.2fjj] . we obtain the desired bound 0(/3^) e"^. This completes the proof of ()3.15|1 . □ 

Proof of fl3.17j) . Since ui H uj2 7^ 0, there is a sequence of bonds bi, . . . ,bn such that ui and uj2 
meet for the first time at b^, share 61, . . . , 6„, and split at 6„ [uJi and 002 may share a bond again 
after 6„). This means that, together with q{u) q{v) q{x — y) q{x — z) in ()3.17|) . the left-hand side of 



15 



()3.17p is bounded by the convolution of two non-vanishing bubbles and YYi=i Qi'^i)'^ ^ where each 
Wi is the spatial component of bi — b^. Using ()3.6p and < ||g||oo, we can bound ()3.17|) by 

e y -^^^^I^ [(1 - ^) V {e\\D\\^)] , < 0{(3') e\ (3.28) 

{ly sy/^ ) \\\ \\oo)\ [1 V (t - s')]^/2 - ^ ' V I 

e<s<s'<t-e 

where, as before, is used up for the Riemann sum approximation. The above estimate can be 
improved to 0(/5^) for oriented percolation, using p.8|) . This completes the proof of ()3.17|) . □ 

3.2 Asymptotics of pl" 

In this section, we compute the asymptotics of the critical point for (unoriented) percolation. We 
follow the strategy in Section 13.11 as closely as possible. However, there are a number of changes 
due to the fact that we have less control of the lace expansion coefficients. For example, the bounds 
on the derivative of Hp with respect to p are not available in the literature, even though in the 
unpublished manuscript [13^, this derivative is computed. To make this paper self-contained, we 
avoid the use of the derivative, which causes changes in the proof. 
We start with some notation. Let 

Tp = snpipD * r;3)(x), T; = sup rf (x). (3.29) 



d 



We will use the following bounds: 

Lemma 3.1. Fix d > 6. For L sufficiently large, and all p < pc, 

Tp < C/3, < 1 + C(3. (3.30) 

We will defer the proof of Lemma 13.11 to the end of this section. 

To compute the asymptotics of tip, we use ()3.3|) and the bound (cf., Proposition 4.1]) 

Til^^ < T;(2TpT;)^^^ (3.31) 
Note that Lemma f3. II together with ()3.5|) and 1)3.311) immediately imply 

p, = l + 0{P). (3.32) 

We now start the proof to improve ()3.32|) one term further. Together with Lemma 13.11 1)3.311) 

proves that the contribution to Y1'n=2^pT ^(/^^)- Thus, we are left to compute n^^ and n^K 
The goal of this section is to prove 



^ oo oo 

= 2 - + OiP'), = D*\o) + nD*-{o) + 0{^% (3.33) 

n=3 n=3 

Using ()3.5|) and ()3.33|) . we arrive at ()1.14j) . Thus, we are left to prove ()3.33p . 
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We again investigate 7r^°^ and fr^^ separately. First, we compute 7r^\ For percolation, we denote 



(0), 



X) 



by {w <^==^ x} the event that w is doubly connected to x. By definition [TT], vr^"^ = X]xgz<^ 
where 

vrf (x)=Pp(o^x)-5o,,. (3.34) 

We wish to use Russo's formula (see, e.g., |8j) to prove that 7r^°' = -itf^ + 0(/5^). However, Russo's 
formula is restricted to events that only depend on a finite number of bonds, so that we will first 
show that Russo's formula may be applied to n^\x). 

Let i?£ = {x G Z"^ : |x| < i}. We note that, since tt^"^ is finite for any p < p^, there is an r < oo 
such that Xli'^B,. '^p K^) = 0(/5^) for any p < pc- In fact, using the BK inequality (see, e.g., 
and the bound Tp^(a;) < for x ^ o ^10, Proposition 2.2]\ we have 

5^ <'(x) < 5^ rp(x)2 < ^vX^f ^ c5^£(<i-i)+2{2-d) _ oir^-") = 0(/?2), (3.35) 

X^Br XI^Br X^Br i>r 

where we assume r = 0{L'^^^^'^~'^^). Let {E in 5^} be the set of bond configurations whose 
restriction on bonds {u,v} with -u, f G 5^. are in Similarly to (jSSg), if i? = 0(L2^/(<^-6)), then 
for any p < Pc we have^ 

^Pp({o ^ x} \ {o ^ x in < 0(/52). (3.36) 

X 

By ()3.35|) - ()3.3fi|) and the mean-value theorem, there is a p G (l,Pc) such that 
= E<c(^) + ^(/^')= E nXo^x in B,) + 0{/3') 

x£Br xeBr\{o} 

= J2 Pi(o^ X in 5^) + (Pe- 1) Yl dpFp{o ^ X in B^) + OiP"^) 

xeBr\{o} xeBr\{o} 

= 7rr + (p,-l) Yl dMo^xinB,) + 0{(3'). (3.37) 

xe_Br\{o} 

We will later identify 7rj°\ and first show that the second term is 0(/?^). Since the event {o <^=^ a; 
in Bfi} depends only on finitely many bonds, we are now allowed to apply Russo's formula to 
obtain 

Y^ dpPp{o <^==^ X in Bn) = D(v — u) Pp((M, v) pivotal for {o •i=?- x in 5^}), 

x£Br\{o} X&Br\{o} {u,v) 

(3.38) 



^In (23 Proposition 2.2], K is of order 0(L~^+'^) with an arbitrarily small number e > 0, and thus is small when 
L is large. Here, we do not care about the dependence of K on L, and will take K = 0{1). 

^The event {o <^=> x} \ {o ^^=> x in 5^} implies the existence of y ^ such that o < — > x, x < — > y and 
y < — > o occur disjointly. Therefore, by the BK inequality and Propositions 1.7(i) and 2.2 of ^U], we obtain 

^ Pp({o ^x}\{o^xin B,}) < ^ Tpix) T,{y - x) T,{y) < c \y\^^~''^+^'-'^ = 0{R'-^). 
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where the factor D{v — u) arises from the derivative of the bond occupation probabihty of {u,v} 
with respect to p, and where a bond is pivotal for o <^==^ x when o <^=^ x in the (possibly modified) 
configuration where the bond is made occupied, while o <^=^ x does not occur in the (possibly 
modified) configuration where the bond is made occupied. 

Since Pc = 1 + 0(/5), and since, by the BK inequality, (j3.38j) is bounded by 

D{v — u) Pp({o < — > u] o \y < — > x} o {o i — ¥ x}) 

x,{u,v) 

< ^ D{v - u) Tp{u) Tp{x - v) Tpix) < p'^Tp < Tp, (3.39) 

x,{u,v) 

so that the second term in ()3.37|1 is 0(/3^). We are left to analyse the first term TifK We follow 
the strategy around ()3.12p . but the details change somewhat. 

Let Sx denote all self-avoiding paths from o to x, and order the elements in Sx in an arbitrary 
way. Then we can write 

7rf = Fi{LUi,iU2 occupied; E^{uji) f] E^{uj2;uji)) , (3.40) 

Li;ina;2=0 

where -E'_<(ti;i) and E-^{uj2]uji) were defined between ()3.11|1 and ()3.12|1 . In words, the event E^{ijJi) 
holds when uji is the lowest occupied self-avoiding walk path from o to x such that there is an 
occupied bond disjoint path from o to x. The event -^^(^2; cui) holds when uj2 is the highest 
occupied self-avoiding walk path from o to x that is bond disjoint from uji. Since Pi is a product 
measure, if we ignore E^{(jJi) fl Ey{uj2; uJi), we obtain 

Pi (cji, ct;2 occupied) = Pi (cui occupied) Pi (c<j2 occupied). (3-41) 

u)inaj2=0 ujinLiJ2=0 

We can then follow the rest of the argument between ()3.14|) and ()3.19|) to arrive at 

*i°' = ^E E W^iM W^iM + 0(/32), (3.42) 

XT^o ijJ\,bJ2&Sx 
|a;i| + laJ2|>3 

where we recall the definition of WpiiS) in (jl.lj) . Here the factor 1/2 has the same origin as the 
one in (|3.18p . and the restriction that \ijJ\ \ + |to'2| > 3 is due to the fact that the smallest cycle 
in percolation has length 3. In ()3.42|) . each iOj is a self-avoiding path from o to x. However, 
as estimated in Section |21 the contribution in which lox or 102 has a self-intersection is 0(/5^). 
Therefore, we can remove the self-avoidance constraint in ()3.42|) . Performing the sum over x ^ o 
and writing uj = {ui, UJ2), which is a random walk path starting and ending at o of length at least 
3, we obtain 

00 

= 2 E (1^1 - 1) + oiP') = 2 - 1) ^*"(«) + (3-43) 
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where \uj\ — 1 = n — 1 is the number of vertices along u, excluding the starting and ending point 
of u. This completes the computation of the leading asymptotics of n^K 

We next derive the asymptotics of tt^^ , following the strategy in [THl IT^ , where the first three 
coefficients of the asymptotic expansion into powers of {2d)^^ of the critical value Pc for nearest- 
neighbour percolation were computed. The details of the argument are changed considerably 
compared to IT^ . Indeed, since we are only interested in the leading order term, while in JH] 
the first three coefficients are computed, many terms that need explicit computation in ^Hl IH] will 
be error terms for us. On the other hand, since in fHl 1^ the asymptotics in nearest-neighbour 
models for large dimensions are considered, long loops lead to error term in IT^ , whereas they 
contribute to the leading asymptotics here. We follow the proof in [THl Section 4.2] as closely and 
as long as possible, and indicate where the argument diverges. 

To define if!^\ we need the following definitions. Given a bond configuration and A C Z*^, we 

recall that x and y are connected through A, and write x < — > y, if every occupied path connecting 
X to y has at least one bond with an endpoint in A. As defined below (I3.7|l . the directed bond 
{u, v) is said to be pivotal for x < — > y, if x < — > u and v < — > y occur, and if x < — > y occurs only 
when {u,v} is set occupied. (Note that there is a distinction between the events {{u,v) is pivotal 
for X < > y} and {{v, u) is pivotal for x < > y} = {{u, v) is pivotal for y < > x}.) Let 

E'{v, X] {u', v') occupied & pivotal for v < — > x s.t. v < — > n'}. (3.44) 

Then, by definition [TT] . 

' = 5^ 5Z - ^) ^« [Mo^u} P, {E'{v, X- &r\o)))\ , (3.45) 

X {u,v) 

where the sum over (m, v) is a sum over directed bonds. On the right-hand side, we use subscripts to 
identify the different expectations. Thus, the subscripts do not refer to the percolation parameter 
p. The cluster C'^"''''(o) appearing on the right hand side of ()3.45|) is random with respect to the 
expectation Eq, but C^"''''(o) should be regarded as a fixed set inside the probability Pi. The latter 
introduces a second percolation model which depends on the original percolation model via the set 
C''^''"\o). We refer to the bond configuration corresponding to the j*'^-expectation as the "level-j" 
configuration. 
By (ESU), 

< tt;^' < 2T;TpT;. (3.46) 

We will use refinements of this bound in the following. 

We first claim that the contribution to ()3.45|1 due to m 7^ o is an error term of order 0(/?^). 
Indeed, if m 7^ o then at level-0 the origin is in a cycle of length at least 3. Standard diagrammatic 
estimates then allow for the replacement in ()3.4(ij) of a factor by a constant multiple of Tp. This 
improves the bound from 0{(3) to 0(/32), by 

We are left with the contribution to ()3.45|) due to -u = o, namely 

Y,PD{V) ¥.,[¥,{E'{v,x-&r\o)))]. (3.47) 
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If X ^ C^°''''(o), then to obtain a non-zero contribution to Fi{E'{v,x;C^°'''\o))), x must be in 
an occupied cycle of length at least 3, in level-1 (in the language of [Si Section 3], the sausage 
containing x must consist of a cycle containing both x and an endpoint of the last pivotal bond 
for the connection from o to x). In this case, in fl3.46|) . we may again replace a factor by a 
constant multiple of Tp, and again this contribution is 0(/3^). We are left to consider 

J^PDiv) Eo[l{xec(-'(o)} P,(^'(i;,a;;C'<°-'(o)))]. (3.48) 

x,v 

This is as far as the analogy with the argument in fSl Section 4.2] goes. We now need to adapt 
the proof there to compute the asymptotics of tt^^ when L — oo. 

If a; G Cf^"'"-' (o) , and if v < — > x, then v ^ > x. We next claim that the intersection with 

the second event in ()3.44|) leads to an error term. We write 



Ux e c^r\o)] iAE'{v,x-c^r\o))] 

= I,[x e C^r\o)] h[v^x] 

X 1 1 — Jj [3(m', f ') occupied & pivotal for v < — > x s.t. v ^ > -u'] 



where we write Jq and /i for the indicator functions on levels and 1, respectively. The latter term 
can be bounded by 

'Ylolze (^^"'(o)] h [{v < — > z}o{z< — > u'} o {(«', v') occupied} o {v' < — > x}] , (3.49) 

(u',v') z 

which, using the BK inequality, yields a bound of the form 

pD{v) Po(o < — > X, o < — > z) Tp{z — v) Tp{u' — z) pD{v' — u') Tp{x — v'). (3.50) 

3:,v,z {u',v') 

By the tree-graph inequality PP 

Po(o < — > X, o < — > <Y Tp{y) Tp{x - y) Tp{z - y), (3.51) 

y 

so that we end up with 

J^ipD* rp)iy) Tpix - y) Tp{z - y) Tp{z -v)irp*pD* rp)(x - z) < = Oi(3^), (3.52) 

x,z,y 

which indeed is an error term. Thus, using the identity 

{x G C'^°''''(o)} = {o < — > X without using (o,t;)}, (3.53) 

we end up with 

^p' = Ep^(^) ^p''''^^) - ^) + ^('^')' (3.54) 



20 



where 

T^°'''\x) = P(o < — > X without using (o, v)). (3.55) 

Note that we can think of t!^°'''\x) as the two-point function on 7/-, where the bond (o, f) is 

removed. We will denote the resulting graph with vertex set 7/ and edge set {{x, ?/} : x,y G 
Z"^, {x,y} 7^ {o, f}} by Z^^^^, so that t^°-''\x) is the two-point function on Z^^^^. We will use this 
observation to compute t^°'''\x). 

We investigate the main term in the right-hand side of ()3.54j) further. Russo's formula, together 
with the BK inequality, yields that 

dpTp{x) = ^ D{z — y) F{{y, z) pivotal for o < > x) < {rp* D * Tp){x), (3.56) 

dpT^"'^\x) = Y,D{z- y) P((y, z) pivotal for o « x in Zf^ ^,) < (r^ * D * Tp){x). (3.57) 
Therefore, we obtain that for p = Pc, 



7r< 



!^ = Y,pMvyp:\x)TpXx -v) + o{(3') 

= J2 D{vyr\x)n{x -v) + 0{{p, - l)TpJ + 0{[3') 

= J2Div)T[°^''\x)n{x-v) + 0{/3^), (3.58) 

x,v 

since Pc = 1 + 0{(3). Furthermore, an argument similar to the one for tt^"' shows that 

Ti{x)=G{x) + 0[{G*g*G){x)), (3.59) 
r{°-'(x) = G{x) (1 - 5,,,) + 5.,.(D*' * G){x) + 0{{G*g* G){x)) +0{D*{G- 6o){x)). (3.60) 

where we recall G{x) = Yl'^=o ^*"' i^) and define 

g{x) = G{x){D*G){x). (3.61) 

We will prove ()3.59p - ()3.60|) in full detail below, and first complete the proof subject to ()3.59|) - 
()3.60|) . Using ()3.59|) - ()3.60|) . together with the fact that for u ^ o, we have G{u) = {D * G){u), we 
end up with 

X^V X 

oo oo oo 

= ^(n - l)D™(o) + ^ Z}*"(o) + 0{f3^) = D*\o) + ^ riZ}™(o) + 0(/3='), (3.62) 

n=2 ?i=3 n=3 

where we use 

{D * 6-*^ * (^)(o) < \\D * G'*3||^ ll^^lli < 0(/32), (3.63) 

for d > 6, by (irTn|) . 

This completes the proof subject to ()3.59|) - ()3.60|) and Lemma IXTl □ 
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Proof of fl3.59p - ()3.60p . We start by proving ()3.59|) . and then adapt the argument to prove ()3.60|) . 
To see fl3.59|) . we recall the arbitrary ordering of the elements in introduced above ()3.40|) . Then 
we have that 

ri(x) = Fp{uj occupied; F^{uj)), (3.64) 

where F^{uj) is the event that uj is the lowest occupied path in Sx- Thus, we can write 

Ti{x) = Pp(co' occupied) — Fp{uj occupied; ^^(u;)'^). (3.65) 

The former term equals 



So,. + (1 - So,,) J2 n + 1) - ^('))- (3.66) 

Clearly, by using inclusion-exclusion on the fact that u is self-avoiding, as in ()2.5|) . ()3.66p equals 

G{x) +0{{G*G){x){G{o) -1)), (3.67) 

which is a contribution to the error in ()3.59|) when we note that G{o) — 1 = [D * G){o). Similarly, 
the second term in ()3.65|) is bounded by 0((G' * g * G){x)^ using the fact that there must exist a 
■u G Z'^ such that there exist bond disjoint occupied paths from a to u, two occupied paths from 
■u to f (of which at least one is non-vanishing) and one from v to x. Thus, by the BK inequality, 
this term is bounded by 

G{u) G{v -u){D* G){v -u)G{x-v) = {G*g* G){x). 

The proof of (|3.6(J|) follows the same ideas. In (|3.65|) and (j3.66|) . we only need to sum over 
self-avoiding walk paths that do not use the bond (o, f). When x = v, this means that |ci;| > 2, so 
that we obtain 

r{°''"(t;) = (D*2 ^ _^ o{{G * g * G){v)) + 0{D{v){G{o) - 1)). (3.68) 
When X ^ v,we can use inclusion-exclusion on the fact that the bond (o, v) is not used, and obtain 

t[''-'\x) = n{x) + 0{D{v)G{x - v)), (3.69) 
and then use (j3.59|) . □ 

Proof of Lemma \3. 11 We use [TT| (5.20)], which states that uniformly in p G [1, Pc) and for L large 
enough 

fp{k) < ^i^, (3.70) 
' ~ 1-D{k) 
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where o(l) tends to when L oo. We also use the standard bound (see e.g. p]) that for x 7^ 0, 

Tp{x) <{pD*Tp){x). (3.71) 

We then follow the proof as in For Tp, we fix x and extract the term in ()3.29|) due to the case 
where every argument of Tp is o, which is pD{x) < pCP (see (|1.8|) ). This gives 

Tp{x)<pC(3 + p Yl Tp{y)Tp{z-y)D{u)Tp{x + z-u). (3.72) 
Therefore, by ()3.71|) . 

Tp<pC(3 + 3p^ sup(D*^ * r;^)(x), (3.73) 

X 

where the factor 3 comes from the 3 factors Tp whose argument can differ from o. In terms of the 
Fourier transform, this gives 

Tp<pCP + 3p' snp [ ^D{k)HpikYe-^'-=pCP + 3p' [ ^ b{kf fp{k)\ (3.74) 

where = [— vr, vr]"^ and we use Tp{k) > [P. For L ^ 1, by ()3.70|1 we obtain 

T,,,C,.v/^^g,^. (3.75) 

Let Bi/i = {k E Dtt : cL^|A;p < r/}. Using ()1.8|1 . we estimate the contribution to the integral in 
dUSD from keU^\ B,/L by 

On the other hand, the contribution from k ^ Bin is, again using ()1.8|) . bounded by 



This proves the bound on Tp. 

The bound on is a consequence of < l + 3Tp. Here the term 1 is due to the contribution to 
()3.29p where the arguments of the three factors of Tp in Tp in ()3.29|) are equal to o. If at least one of 
these arguments is nonzero, then we can use ()3.71|) for the corresponding two-point function. □ 

A Bounds on D*'^{x) and q^'^^x) 

In this appendix, we prove ()3.6|) for any e G (0, 1] assuming p.Sj) . The inequality p.lO|) follows 
by taking e = 1. 
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First, we note that 

n-l 



= (1 - er6o,. + 5^(1 - * q*^){x), (A.l) 

j=0 

where we suppose that the empty sum equals zero. When n < N = e^^, we use (IA.1|) to obtain 



n-l 



0(/3) 



< (1 - er6o,. + 5^(1 - eYe \\D\\^ < {1 - sTdo,, + 0{(3) < (1 - + , 

^ [IV 

(A.2) 

as required. On the other hand, when n > N, we use 

g*"(a;) = (1 - e)"^,,, + (1 - + ^^(n - 1 - j)(l - £)"-2-.^2(^*2 ^ (^ g) 

i=o 

which is obtained by substituting ()A.1|) into g*-' in the right-hand side of ()A.1|) . Since the second 
term in the right-hand side is bounded by 0{P)nee~^^ < 0{(3) [lV(ne)]~'^/^, it suffices to investigate 
the third term. Let 5*1 be the sum over j < N, and let 5*2 be the remaining sum, i.e., 

S,= in-l-j){l-er-'-^e\D*Uq*n{x), (A.4) 

0<j<N 

52= J2 {n-l-j){l-e)''~^''e\D*Uq*^){x). (A.5) 

N<j<n-2 

For Si, we use ()A.2|) to obtain 

n— 1 n— 1 

^1 < (1-£)""VD*'(o) ^ + 0{f3)e^ J2 - ^)'"' < ^(/S) (ne) V"^ (A.6) 

l=n-N l=n-N 

For 5*2, we recall the definition Bi/l = {/c G Dtt : cL^|/cp < t]} below ()3.75p . and let B^^^ = {k E 
Bi/L ■■ qik) > 0}. Note that 

("•^ <- 4^ S ^ Lt,. Sr ^"='^""='1 ■ 

Recall (jl.Sp . For the first integral, we use 

g(A;) = 1 - £[1 - D{k)] < e-^ti-^^'^)! < e-^^^'I'^l', (A.8) 

while for the second integral in ()A.7|) . we use, noting that without loss of generality, we may assume 
that r] < 1, 

\q{k)\ <\l-e{2-ri)\y (l-eri) <l-6ri. (A.9) 
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Therefore, we have 

= 0(/5) + 0(/3) (1 - er^y. (A.IO) 

Substituting ()A.10|) into ()A.5|) and separating the sum of the first term in ()A.10|) into the sum 
over N < j < — 1 and the sum over f — 1 < j < n — 2, we obtain 

S,<0{(3){ne-e-l){l-et^-^e ^ {je)-"!^ 

Af<j<f-1 

+ 0(/3) (f - 5^ ^ _ 1 _ j) (1 _ 

t-l<i<n-2 

+ 0(/?)(l-er/r-2e^ 5^ (ri-l-j) 

< 0{(3) (n£)e-"^/' + 0(/3) {ne)-''/^ + 0(/3) (n^)^' e^''^^ (A.ll) 

The proof of (jSIHl) is completed by combining (E^-dQ, (|^^ and (|XTT| . and (n£)2e-"^/2 < 
C[l V for ^ > □ 



B Computation for the spread-out uniform model 

In this appendix, we compute the model-dependent terms of — 1 in p.l2|) - p.l4|) when the 
probability distribution D is defined as in Recalling p.l5|l . we have 

D*\o) = ^^^_^\^,_^ = I + 0{f3L-') = (3 W\o) + 0{(3L-'). (B.l) 

This relation can be extended as follows: 

Proposition B.l. Let D be the function defined in (|1.9|) . For a = 0,1, as L ^ oo, 

oo oo 

J2in + 1)" D*-{o) =PY.^n+ ir W-{o) + 0(/?L-^), (B.2) 

n=3 n=3 

oo oo 

^Z}*2"(o) =/3 5^t/^2"(o) + 0(/3L-i), (B.3) 

n=2 n=2 

where d> A + 2a m and d > 4 m (P3|) . 

Theorem 11.21 is an immediate consequence of ()B.1|) and Proposition IB. II Note further that the 
coefficients of (3 in ()B.2|) - ()B.3|) are bounded if d > 2 + 2a; and d > 2, respectively, which suggests 
that also for = 3 + 2a and = 4 + 2a the leading order contributions should be given by the 
first terms in (lR2l)-(lR3l). 
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Proof. For x G Z'^, define 



Do{x) 



l{\\x\\oc<L} 

{2L + iy' 



(B.4) 



to be a regularized version of D in ()1.9p . It is obvious tliat D*'^{o) satisfies the same estimate as 
in dnUl). We note that 



a:::0<||a;||oo<-Z' 



(2L + 1)'^[(2L + l)-^ - 1] 



{2L + IY 



^^^^ [(Z^*^- * D<--^^){o) - {D<^-'^ * /^:(™-^'))(o)J . (B.5) 



By this identity and the fact that D*^{x) also satisfies p.lO|) . we can approximate the expressions 
in the left-hand side of ()B.2jl - ()B.3|l by the corresponding expressions defined in terms of Do instead 
of Z), up to 0(/?^) when > 4 + 2a and d > 4, respectively. For example, for ()B.2|) with « = 0, 
we use ()1.10|1 to obtain 



Y.d*-{o)-y,d: 



?i=3 



n=3 



oo oo 



j=l n=jV3 



= J^ZTW ^ '^"^ * + * Do + * G) * Go) (o) = 0(/?2), (B.6) 

where 5o(a;) = So,x and Ga(x) = J2'^=o i^) ■ Therefore, to prove Proposition IB.lf it suffices to 
show that, for a = 0, 1, 



oo oo 



n=3 



(B.7) 
(B.8) 



n=2 



71=2 



We prove ()B.7|) for a = by comparing the Fourier transform of the first term in the right-hand 
side of (jRTjl . i.e.. 



d'^k U{kf 



-U{k) 



(B.9) 
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with the Fourier transform of the left-hand side of ()B.7|) . i.e., 



n=3 



where Df 



and also 



and 



(B.IO) 



U{k) 



d'^x U{x) e 



ik-x 



nsin kj 



{2L + 1) sin f 



(B.ll) 

(B.12) 
(B.13) 



The simple product formula ()B.13|) is the main advantage of using Do instead of Z^. It follows 
from ()B.13|) that Do{-[^) = U{k)/U{^^^), which approximates U{k) for large L. We write 



(B.14) 



n=3 



n=3 



n=3 



where, by ()B.11|) . the first term is 0{[3L ^), and 



h 



d'^k U{kf 



d'^k 



1. (2^)' 



U{kf 



k \3 



(B.15) 



We prove below that each Ij is 0{L~^) ii d > 2. This suffices to prove ()B.7|1 for a = 0. If fact, we 
will prove that each Ij is 0{L~'^ logL) if d > 2, which also identifies the coefficient of l3L~^. 

To estimate each Ij, we use the following properties of U{k) and Do[j^) that follow from the 
standard estimates for the trigonometric functions: for any k, 



mk)\<U^iy\k,\)-\ 



l-f>(A:)>Ci(lA|A;n, 



(B.16) 



and for k E D(^^_^_i^^, 



\Do{lTT)\<C2lliiy\k,\r\ 



l-^o(^) >C3(lA|A;p 



\U{k) - Do{j;^) \ < c,L-'\U{k)\\k\' 



(B.17) 
(B.18) 
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where each q G (0, oo) is independent of L and k. To see the first inequahty in ()B.16|) . we only 
need to use |^^| < (1 V r)~^ for any r. For the first inequahty in ()B.17|) . we recaU Z)o(-^-^) = 



U{k)/U{ 



2L+1 



and use sisj: > 1 for any r G [-|, |], so that l^^lf^)"^! < {^Y, and then use the 



bound on \U{k) \ in (jRT6|l . The second inequahties in (fRlfijl - ffRTrjl foUow from (fRT2|l and jEHl), 
respectively. Finally, for ()B.18|) . we again use ^ < < 1 for any r G [— f , f ] to obtain 



it/(fc)-^,(^)i<(l) \uik)\ n 



sm 



2L+1 



sin^ 



j=l 2L+1 j=l 2L+1 

which is bounded by \U{k)\ 0{\k\'^L-^), using < ^ - 1 + J < ^(f )2 for any r G [-f , f]. This 
completes the proof of (lRT6l) - (lBl8l) . 

First, we consider Ji. Using ()B.16j) . we obtain 

d 



\h\ < c 



For I2, we use 



(fk \lJ{k)\^ < c'E 



(L + i),r 



1 L^{L+|)7r kf 



n 



;iv|fc,l)= 



0{L 



-2\ 



rf3 



1 — M 1 — (i 



\u — d\ + ud + (P ~ ud{u + d)\ ^ \u — d\ {v? + ?>\ud\ + 



[l-u){l-d) 



with M = f/(A;) and d = bo{j^). Using (fRJ6|) - (fRT8|l . we obtain 



I/2I < cL 



-2 



d^k 



fepn;=i(ivife.i) 

(1 A |A;|2)2 



-3 



O(L-^logL), 



d'^k ^"(lTx)^ 



n=2 -^^{L+h 



This completes the proof of Proposition IB.l 



(B.20) 



(B.21) 



(B.22) 



for d>2. The proof of (lR7|l for a = is completed by (iBlll . (Ib:20|) and (|R22l). 

The same strategy explained above also applies to the proof of (jB.2jl for a = 1 and ()B.3jl , using 
the following expressions: 

d^k MijiYi^-^Miir)] 

1 2 ' 



(B.23) 

(B.24) 
□ 



Acknowledgements 

The work of AS was supported in part by NSERC of Canada. The work of RvdH and AS was 
supported in part by Netherlands Organisation for Scientific Research (NWO). This project was 
initiated during an extensive visit of RvdH to the University of British Columbia, Vancouver, 
Canada. We thank Gordon Slade for comments on a preliminary version of the paper, and for 
pointing us to the appropriate bound in ^T] that implies Lemma f3. 11 



28 



References 



[1] M. Aizenman and C. M. Newman. Tree graph inequalities and critical behavior in percolation 
models. J. Statist. Phys. 36 (1984): 107-143. 

[2] D. J. Barsky and M. Aizenman. Percolation critical exponents under the triangle condition. 
Ann. Probab. 19 (1991): 1520-1536. 

[3] C. Bczuidcnhout and G. Grimmett. Exponential decay for subcritical contact and percolation 
processes. Ann. Probab. 19 (1991): 984-1009. 

[4] M. Bramson, R. Durrett and G. Swindle. Statistical mechanics of crabgrass. Ann. Probab. 
17 (1989): 444-481. 

[5] C. Borgs, J.T. Chayes, R. van der Hofstad, G. Slade, and J. Spencer. Random subgraphs of 
finite graphs: II. The lace expansion and the triangle condition. Preprint, (2003). 

[6] J. T. Gox and R. Durrett. Oriented percolation in dimensions d > 4: bounds and asymptotic 
formulars. Math. Proc. Cambridge Philos. Soc. 93 (1983): 151-162. 

[7] R. Durrett and E. Perkins. Rescaled contact processes converge to super-Brownian motion in 
two or more dimensions. Probab. Th. Rel. Fields 114 (1999): 309-399. 

[8] G. Grimmett. Percolation. Springer, Berlin (1999). 

[9] G. Grimmett and P. Hiemer. Directed percolation and random walk. In and Out of Equilibrium 
(ed. V. Sidoravicius). Birkhauser (2002): 273-297. 

[10] T. Hara, R. van der Hofstad, and G. Slade. Critical two-point functions and the lace expansion 
for spread-out high-dimensional percolation and related models. Ann. Probab., 31 (2003): 
349-408. 

[11] T. Hara and G. Slade. Mean-field critical behaviour for percolation in high dimensions. 
Commun. Math. Phys., 128 (1990): 333-391. 

[12] T. Hara and G. Slade. The self-avoiding-walk and percolation critical points in high dimen- 
sions. Combin. Probab. Comput., 4 (1995): 197-215. 

[13] R. van der Hofstad. The derivative of the lace expansion coefficients for unoriented percolation. 
Unpublished document (2003). 

[14] R. van der Hofstad and A. Sakai. Gaussian scaling for the critical spread-out contact process 
above the upper critical dimension. Preprint (2003). To appear in Electr. Journ. Probab. 

[15] R. van der Hofstad and G. Slade. A generalised inductive approach to the lace expansion. 
Probab. Th. Rel. Fields 122 (2002): 389-430. 



29 



[16] R. van der Hofstad and G. Sladc. Convergence of critical oriented percolation to super- 
Brownian motion above 4+1 dimensions. Ann. Inst. H. Poincare Probab. Statist. 39 (2003): 
413-485. 

[17] R. van der Hofstad and G. Slade. The lace expansion on a tree with application to networks 
of self-avoiding walks. Adv. Appl. Math. 30 (2003): 471-528. 

[18] R. van der Hofstad and G. Slade. Expansion in for percolation critical values on the n-cube 
and Z": the first three terms. Preprint (2004). To appear in Combin. Probab. Comput. 

[19] R. van der Hofstad and G. Slade. Asymptotic expansions in for percolation critical values 
on the n-cube and Z". Preprint (2004). 

[20] H. Kesten. The critical probability of bond percolation on the square lattice equals |. Com- 
mun. Math. Phys. 74 (1980): 41-59. 

[21] T. Liggett. Stochastic Interacting Systems: Contact, Voter and Exclusion Processes. Springer, 
Berlin (1999). 

[22] T. Liggett. Stochastic models of interacting systems. Ann. Probab. 25 (1997): 1-29. 

[23] N. Madras and G. Slade. The Self-Avoiding Walk. Birkhauser, Boston (1993). 

[24] B. G. Nguyen and W.-S. Yang. Triangle condition for oriented percolation in high dimensions. 
Ann. Probab. 21 (1993): 1809-1844. 

[25] M. D. Penrose. Self-avoiding walks and trees in spread-out lattices. J. Statist. Phys. 77 
(1994): 3-15. 

[26] A. Sakai. Mean-field critical behavior for the contact process. J. Statist. Phys. 104 (2001): 
111-143. 



30 



